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ABSTRACT 


On 24 April 1975, numerous severe storms developed in an area 
extending from southwest Oklahoma to eastern Tennessee during the after- 
noon and evening hours Significant severe weather reports included 
2.75 in. hail at Wewoka, Oklahoma and a destructive, killer tornado at 
Neosho, Missouri An unusually comprehensive data set was available on 
this day that included 5-rainute interval SMS visible and infrared 
radiances from 1800 to 0200 GMT, intensity contoured radar reflec- 
tivities, special upper-air analyses from NASA's Atmospheric Variability 
Experiment, and low-level wind fields from satellite inferred cumulus 
velocities. 

The mesoscale nature of the forcing and evolution of these storms 

IS investigated, with emphasis on techniques to aid in the early 

detection of such severe events. In the pre-storm environment (t-4 to 

t-2 hours), the satellite wind fields were combined with moisture 

parameters to derive horizontal moisture flux information Low level 

“3 -1 -1 

moisture convergence of 1 5 x 10 g kg s or greater was indicative 
of regions of subsequent severe storm genesis. Dynamic parameters 
such as boundary layer vorticity production and relative vorticity were 
also useful prognosticators of subsequent severe activity 

In the first 30 minutes of cloud lifetime, infrared (IR) cloud 
-3 -1 

growth rates of 2 x 10 s or greater accurately indentified potential 
severe storms Little success was realized in forecasting the exact 
time of the severe event on the ground. A comparison of IR determined 
cloud tops with radar observed tops illustrated particular problems 


11 



in inferring cloud height for storms penetrating the tropopause. The 
need for rapid-time interval satellite coverage is stressed, as is the 
need for a computer-linked video system to process and interpret data. 



TABLE OE CONTENTS 


Chapter Page 


LIST OF TABLES Vlll 

LIST OF FIGURES ix 

1.0 INTRODUCTION 1 

2.0 CASE STUDY; 24 APRIL 1975 5 

2.1 Synoptic Situation 5 

2.2 Severe Storm Reports 12 

3.0 MESOSCALE WINDS AND MOISTURE 15 

3.1 A Wind Tracking System - AOIPS 15 

3.2 Satellite Derived Wind Fields 18 

3.3 Moisture Convergence 22 

3.3.1 Moisture convergence using satellite 

derived precipitable water 26 

3.3.2 Moisture convergence from surface mixing 

ratios 30 

3.4 Dynamic Parameters . , , . . ,33 

3.4.1 The divergence term 33 

3.4.2 Relative vorticity 36 

3.5 Summary of Wind and Moisture Fields 40 

3.5.1 Limitations 40 

4.0 SATELLITE MONITORING OF SEVERE LOCAL STORMS 43 

4.1 Previous Studies 43 

4.2 Radiative Transfer Equation 47 

4.3 Data and Methods 48 

4.4 Results for Severe Storms 50 

4.4.1 Wewoka 50 

4.4.2 Neosho 54 

4.4.3 Cloud no. 1 57 

4.5 Comparison with Non--Severe Storms 57 

4.6 Results for All Storms 62 

4.7 Summary of Infrared Growth Rates 65 

5.0 RADAR-SATELLITE COMPARISON 67 

5.1 Areal Comparisons 67 

5.2 Storm Top Comparisons 

6.0 CONCLUSIONS 76 

REFERENCES 79 

ACKNOWLEDGEMENTS 82 


iv 



Chapter 


Page 


APPENDIX A- -ALL DIGITAL VIDEO IMAGING SYSTEM 

FOR ATMOSPHERIC RESEARCH 83 

APPENDIX B— GLOSSARY OF TERMS .87 


V 



LIST OF TABLES 


Table 


Page 

1-1 

Objective Indicators o£ Severe Storm Activity 
(from Endlich and Mancnso, 1968) 

4 

3-1 

Statistics on Cloud Location and Motion . ... 

19 

3-2 

Previous Studies of Severe Weather-Related 
Parameters . . . . ... ... 

41 

4-1 

Convective Rainfall Estimation Scheme (from 
Scofield and Oliver, 1977) . . ... 

44 

A-1 

All Video Imaging System for Atmospheric Research 
(ADVISAR) Hardware Components . . . 

85 


vx 



LIST OF FIGURES 


Fxgure Page 

2-1 Surface analysis 24 April 1975 1600 GMT . . 6 

2-2 Surface analysis 24 April 1975 2100 GMT . . . .7 

2-3 500 millibar analysis 24 April 1975 2100 GMT . 9 

2-4 300 millibar analysis 24 April 1975 2100 GMT . 10 

2-5 Monett, Missouri upper air sounding 24 April 1975 

2100 GMT . . . . . 11 

2-6 700 millibar vertical velocity in ybars s ^ (from 

Wilson, 1976) Regions where severe thunderstorms 
developed are cross-hatched .13 

2- 7 Severe storm reports 24-25 April 1975 . . . 14 

3- 1 SMS-2 visible channel image at 2058 GMT . . . 16 

3-2 A portion of the AOIPS used for cloud tracking . 17 

3-3 Satellite derived wind field at 1800 GMT . 20 

3-4 Satellite derived wind field at 2100 GMT 21 

3-5 Objectively analyzed wind field at 1800 GMT . . .23 

3-6 Objectively analyzed wind field at 2100 GMT . . 24 

3-7 (from Sasaki, 1973) a) Surface moisture convergence 

(10"^ g kg'^s"^) 19 April 1972 1800 CST, b) NWS radar 
summary 19 April 1972 1945 CST . 25 

3-8 Satellite derived precipitable water field (cm) 

at 1615 GMT . . ... . . 27 

3-9 The moisture divergence field at 1800 GMT Only 

negative values (moisture convergence) have been 

contoured Units are cm H2O s”^ x lO'^ . . 29 

3-10 NWS radar summary at 1835 GMT ... 29 


3-11 2100 GMT moisture divergence field using satellite 

derived winds Only negative values (moisture 
convergence) have been contoured Units are 
-1 -1 -5 

g kg s X 10 and contour interval is 

50 X 10 ^ g kg ^s ^ • . 31 


vii 



Figure 


Page 


3-12 

3-13 


3-14 


3-15 


3-16 


4-1 


4-2 


4-3 


4-4 


4-5 


4-6 


4-7 


Severe storm locations at 2300 GMT .... 32 

2100 GMT moisture divergence field using observed 
surface winds Units and contour interval are the 
same as in Figure 3-11 . 34 

"Vorticity production" at 2100 GMT Only positive 

-2 -9 

values have been contoured Units are s x 10 

-9 -2 

and contour interval is 5 x 10 s .... .37 


Relative vorticity computed from satellite derived 
winds at 2100 GMT Only positive values have been 

-1 -4 

contoured Units are s x 10 and contour 

interval is 5 x 10 ^s ^ . . . 38 


Relative vorticity computed from the modified 

satellite derived winds. Details are as in 

Figure 3-15 ... 39 

Thunderstorm anvil models (from Auvine and Anderson, 

1972) a) A model of airflow within the thunderstorm 
anvil assuming calm environmental winds, b) The 
material outline of a cumulonimbus anvil assuming 
Its growth to be the result of a point source in a 
uniform stream . . 46 

SMS- 2 infrared image at 2242 GMT Clouds examined 
in this study are labelled . . .49 

Thunderstorm growth rate diagram for the Wewoka 
hailstorm. N is the number of data points in the 
element with blackbody temperature T^o < T , where 

DD 1 

T^ IS as identified in the diagram .51 

Temporal change in divergence and rate of decrease 

of cloud top minimum temperature (top) for the 

Wewoka storm . . ... ... ... 53 


Thunderstorm growth rate diagram for the Neosho 
tomadic storm . 55 

Temporal change in divergence and rate of decrease 

of cloud top minimum temperature (top) for the 

Neosho storm . .... 56 

Thunderstorm growth rate diagram for cloud ffl, a 

hail producer 58 


viii 



Figure 


Page 


4-8 

4-9 

4-10 

4- 11 

5- 1 

5-2 

5-3 

5-4 

5-5 

A-1 


Thunderstorm growth rate diagram for several 

non -severe storms . .... 60 

Temporal variation of minimum blackbody temperature 
(Tgg) for two severe and two non-severe storms . 61 

Relation of cold area rate of growth to cloud top 
relative temperature and the occurrence of severe 


weather (from Adler and Fenn, 1977) . . 63 

Results from this study using only the first 
thirty minutes of cloud lifetime . .... 

64 

Intensity contoured radar reflectivities for the 

Wewoka storm . 68 

Character display of the infrared growth rate of 

the Wewoka storm ... .69 

Radar echo growth rates for the Wewoka storm . 70 

Comparison of radar and satellite indicated 

cloud top heights .... 73 

Conceptual model of satellite and radar 

viewed cloud area 75 

Ml ^igital ^ideo imaging ^stera for Atmospheric 
Research (ADVISAR) 85 


IX 



1.0 INTRODUCTION 


With the laimch of geostationary satellites in the late 1960 *s, 
a new era was introduced in the prediction, detection, and understand- 
ing of the atmospheric phenomena known as severe local storms. A 
sequence of satellite images provides one with an immediate feel for 
the dynamics of the situation one is observing. The combination of 
satellite measurements with conventional surface and upper air 
observations provides a more complete description of "weather" than 
either method alone 

It IS the piurpose of this study to use such a combination of data 
sources in the development of techniques to better predict the occur- 
rence and severity of severe local storms. A working definition of 
the phenomena is first needed 

Revere I^ocal ^torms (SELS) dangerous storms that usually cover 
relatively small geographic areas or move in narrow paths and are of 
a sufficient intensity to threaten life and property.* 

For the purposes of this study, severe events included tornadoes, 
funnel clouds, and thunderstorms with either hailstones greater than 
0 75", or wind gusts greater than 50 knots The aim of the prediction 
scheme to be described is to provide a "point" forecast one to two 
hours prior to the onset of a severe storm. The technique requires 
high temporal frequency satellite data (5-15 minutes) and a computer- 
linked video system that allows for man-machine interaction Both 
were available in this "after-the-fact" diagnostic study. 


*National Severe Local Storms Operations Plan, U.S. Department of 
Commerce — NOAA, Washington, D.C., 1977, 38 p. 
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The day chosen for study was 24 April 1975, when several severe 
storms developed across an area from southwestern Oklahoma to eastern 
Tennessee during the afternoon and evening hours. Significant reports 
included 2,75'* hail near Wewoka, Oklahoma and a destructive tornado 
at Neosho, Missouri As will be demonstrated in later sections, the 
outbreak of these large convective storms was due in large part to 
mesoscale (meso-3), not synoptic scale forcing. The infamous "super- 
outbreak” of tornadoes on 3-4 April 1974 is a classic example of a 
large mesoscale C^eso-a) development controlling the outbreak of these 
storms. Among the recommendations of a panel on severe storms formed 
soon after this outbreak were 

1^ better detection and tracking of convective activity 

2) development of better statistical relationships between 
convective storms and the large mesoscale disturbances 
that produce them 

3) inprovement in the dynamical prediction of these large 
mesoscale disturbances* 

t 

The first and third recommendations are addressed directly in this 
study 

The monitoring of the pre-storm environment is of utmost inportance 
in determining areas of subsequent severe weather. One type of 
monitoring can be accomplished through the use of low-level wind 
fields derived from satellite determined cumulus velocities. ' These 
fields can provide horizontal moisture flux information when coupled 
with an appropriate moisture parameter (mixing ratios, precipitable 

*Severe Storms-Prediction, Detection and Warning, National Academy 
of Sciences, Washington, DC, 1977, 78 p. 
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water) . ' Storms must necessarily grow in areas of moisture convergence, 
and one attempt of this study will be to compare the magnitude of such 
convergence with subsequent storm severity The dynamic parameters of 
boundary layer vorticity production and relative vorticity will .also 
be investigated with respect to severe storm activity Table 1-1 is a 
comprehensive list of objective indicators of severe storm activity. 

Continued remote sensing of the cloud field in its growth phase 
IS then accomplished through the use of satellite infrared radiances. 
Rates of growth for several storms of differing severity are examined 
in an attempt to develop criteria for a short-term predict lon-waming 
scheme. The problem becomes challenging because direct identification 
of tornadoes or regions of heavy precipitation is not possible when 
viewing the storm from above, the internal structure of the storm is 
hidden from both optical and infrared sensors. Fujita (1973) has 
shown that the height of cloud- tops can undergo fluctuations indicative 
of the changes in the intensity of updrafts within the storm. It is 
this tjrpe of indirect approach that will be undertaken here. 

A final section on the comparison of satellite and radar data for 
severe storms will detail some particular problems in infrared remote 
sensing of cloud dynamics. Limitations in the detection of small, 
intense echo areas and overshooting domes are discussed, as is the 
investigation of radar and satellite determined cloud heights for 
storms penetrating the tropopause. 
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Table Objective Indicators of Severe Storm Activity 

(from Endlich and Mancuso, 1968) 


Quantity 

Rating 

Remarks 

High wind speed in low levels. . . 

Poor. . 

These conditions are 

High temperature and humidity 
in low levels 

Poor . . 

generally present in 
severe storm areas. 

Large moisture gradient near 
700 mb 

Poor . . 

but also cover exten- 
sive regions without 

High wind speed in upper 
troposphere 

Poor . . 

severe storms. Thus, 
they appear inferior 

Large temperature gradient at 
tropopause level 

Poor . . 

to other quantities as 
indicators . 


Unstable lapse rate between low 


and midtroposphere 

Unstable Showalter index 

Fair . . 
Fair . . 

The most unstable lapse 
rate and most unstable 
index lie generally 
south of severe storm 
areas • 

Area between low-level temperature 
and moisture axes 

Fair. . 

These axes are not 
always well developed. 


Destabilizing temperature and 
advection between low and mid- 
troposphere * 

Good * * 

Proximity to 4260 m - thickness 
line (850-500 mb) and to thick- 
ness ridge. . 

Good. . 

Area of negative vorticity of wind 
shear vector between low and 
midtroposphere 

Good . . 

Horizontal divergence in upper 
troposphere 

Good . . 

Upward vertical motion at top of 
boundary layer 

Good. . 

Negative values of V*MV and of 
V*TV in low levels 

Good. . ^ 

Vorticity production in low 
levels 

Good . . ^ 

Fr ontogenesis of temperature and 
of moisture in low levels .... 

Good . . 

Destabilizing distribution of 
V*TV between low and mid- 
troposphere 

Good . . 


^Examined in this study 


2 0 CASE STUDY 24 April 1975 
2 1 Synoptic Situation 

On 24 April 1975, numerous severe storms occurred in a region 
stretching from southvvest Oklahoma to eastern Tennessee during the 
afternoon and evening hours Figure 2-1 is the surface analysis at 
1600 GMT (1000 LST) A cold frontal boundary in the Oklahoma Panhandle 
extended northeastward through Kansas and Missouri,, where it became 
stationary An intense moisture gradient existed across a dry line in 
western Oklahoma. Extremely dry air (dew points 28°-48°F) was moving 
slowly westward from the Texas plains Horizontal wind shear across 
this dry line enhanced the convergence of moist air (dew points dS"- 
70®F) in' eastern Texas and Oklahoma. A dissipating squall line from 
the previous day's convective activity separated hot, moist air in the 
warm sector from cooler air produced by nocturnal thunderstorm 
activity. 

The evolution of the surface system is portrayed in Figure 2-2, 
the surface analysis five hours later at 2100 GMT. The western portion 
of the cold front had moved a considerable distance southward under the 
influence of 20-25 kt northerly winds. Light winds in the air mass 
behind the dry line caused only a slight eastward shift of that boundary. 
South of the front, a narrow tongue of hot, dry air had advanced north- 
eastward into a low pressure area forming on the Oklahoma- Kansas border 
This type of sub-synoptic scale low, shown in Figure 2-2, frequently 
occurs with tornado -producing storms (Tegtraeier, 1974) , and this case 
was no exception. The nocturnal thunderstorm boundary is still evident 
in the temperature field and is shown as a dissipating warm front 
Boundaries of this type have been shown to be important in severe storm 
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Figure 2-1 Surface analysis 24 April 1975 1600 GMT. 
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genesis when they intersect frontal boundaries (Purdom, 1976a) In the 
warm sector, temperatures had climbed into the upper 80 's with dew point 
temperatures in the high 60 's_ 

A comprehensive set of upper air data was available every three 
hours this day from NASA's fourth Atmospheric Variability Experiment 
(AVE IV). The 500 mb analysis at 2100 GMT (Figure 2-3) showed a weak 
shortwave trough extending south from Kansas into central Texas The 
upper air flow field split into two regimes west of the Rocky Mountains. 
A stronger shortwave, within the northern stream, finally merged with 
the southern wave by 0000 GMT on the 25th Although severe storm 
activity occurred throughout the night, the most intense storms were 
those associated with the southernmost wave during the period 2300- 
0200 GMT on the 24th through 25th 

The 300 millibar analysis at 2100 GMT is portrayed in Figure 2-4 
A shortwave trough centered in South Dakota and Nebraska was advancing 
southeastward into a region of upper level difluence This difluent 
region was located between a polar jet stream of 70-80 knot winds and 
a subtropical jet of 60-70 knot winds. VVhitney (1977) has found that 
this type of difluent configuration ahead of an advancing trough is 
favorable for severe storm development IVhitney found that in the 
April 24 case, as in other cases, the greatest potential for release of 
conditional instability existed north of the subtropical jet. 

The Skew T/Log P plot of the upper air sounding taken at Monett, 
Missouri (UMN) at 2100 GMT is shown in Figure 2-5. The sounding is 
typical of a classic severe storm sounding, with a low level moisture 
source capped by an inversion and dry air aloft In the boundary layer, 
mixing ratios were almost 12 g kg and Veering winds indicated low 
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Fxgure 2-3 500 mxllibar analysxs 24 Aprxl 1975 2100 GMT 
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Figure 2-4 300 millibar analysis 24 April 1975 2100 GMT. 
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Figure 2-5 Monett, Missouri upper air sounding 24 April 1975 2100 GMT. 
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level warm air advection. The sounding's lifted index (LI) was -4, 
indicative of the substantial buoyancy a subcloud parcel of air would 
have if lifted to 500 mb. A large area of negative buoyancy is also 
apparent in the 200 mb layer between the lifted condensation level (LCL) 
and the level of free convection (LFC) . 

Using the AVE IV sounding data, Wilson (1976) computed the synoptic 
scale vertical motion field (to) using a kinematic method. His 700 mb 
to field at 2100 GMT is shown in Figure 2-6, with units of ybars/sec 
Regions where severe storms developed are cross-hatched. The large- 
scale vertical motion fields apparently did not organize the convection 
on this day, storms developed in regions of diagnosed subsidence or of 
only slight upward motion The mesoscale nature of the dominate 
forcing mechanisms is discussed in subsequent chapters 

2 2 Severe Storm Reports 

Due to a combination of synoptic scale features (cold front, 
advancing upper level trough, upper level difluence) and mesoscale 
organization (subsynoptic low, dry line) several severe storms developed 
in the afternoon and evening hours of 24-25 April Reports from Storm 
Data* have been summarized in Figure 2-7. The operational range of 
radar sites at Little Rock (LIT), Oklahoma City (OKC), and Kansas City 
(MKC) , are also shown. Radar data from crucially located Monett, 
Missouri were not recorded. Significant storm reports included 
hailstones of 2.75 in. (7.0 cm) diameter at Wewoka, Okhaloma at 2315, a 
tornado near Blue Jacket, Oklahoma at 0000 GMT; and a tornado at Neosho, 
Missouri at 0040. This later tornado caused three deaths, 22 injuries 
and $10.5 million in property damage. 

*A monthly publication of the Environmental Data Service of the National 
Oceanographic and Atmospheric Administration. 
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Figure 2-7 Severe storm reports 24-25 April 1975- 




3.0 MESOSCALE WINDS AND MOISTURE 

3.1 A Wind Tracking System - AOIPS 

Five minute interval SMS visible channel data of 1 km resolution 
were used to derive low-level wind fields by tracking small cumulus 
clouds. Winds were tracked on NASA's Atmospheric and Oceanographic 
Information Processing System (AOIPS)*. A complete description of the 
AOIPS and cloud tracking techniques is given by Billingsley (1976). 

The day chosen as a case study was 24 April 1975, in an area approxi- 
mately 10 degrees on a side centered near Oklahoma City, Oklahoma. 
Intense convective development occurred 2-3 hours after the tracking 
time. A time series of rapid time interval (5 minute) SMS digital 
data was selected; for example, a sequence comprised of the 
2057-2102-2107-2112 GMT images (see Figure 3-1) . These images were 
navigated (by a landmark registration technique), processed into a 
television loop, and displayed on a Hazeltine image processor (Figure 
3-2) . A specific cloud element was then selected and an electronic 
cursor was used to specify the cloud's location throughout the loop. 

The cursor was positioned near the rear of the cloud (relative to its 
movement) to avoid propagation and shear contamination of the velocity 
field. The velocity of the cloud element is computed from its location 
on successive images, knowing the time interval between images. 

The acceptance of a wind vector into the data set was based on 
two considerations: 


*Colorado State University has an equally sophisticated system, the 
All Digital Video Imaging System for Atmospheric Research (ADVISAR) . 
See Appendix A. 
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Figure 3-1 SMS-2 visible channel image at 2058 GMT. 


ORIGINAL PAGE IS 
OF POOR QTJAT.ITY 





Figure 3-2 A portion of the AOIPS used for cloud tracking. 
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1. On an objective basis, the statistics on cloud motion 

(Table 3-1) are examined for consistency and representativeness, e.g., a 
cloud element should not be accelerating greatly. 

2. On a more subjective basis, the vector is shown superimposed 
on the loop, again in an attempt to determine its representativeness of 
the cloud motion. 

A "complete” data set consisted of 100-300 vectors. 

3.2 Satellite Derived Wind Fields 

Two complete wind sets tracked on the AOIPS are shown graphically 
in Figure 3-3 (1800 GMT) and Figure 3-4 (2100 GMT). Lengths are scaled 
such that a vector of length one degree represents a wind speed of 
17.8 ms ^ u velocity and 22.2 ms~^ v velocity. Voids exist in the data 
where either high cirrus clouds or mid-level altocumulus obscured the 
low-level cumuli, or in the dry stable air where convection was 
inhibited. In Figure 3-4, the frontal region running E-W through Kansas 
and Missouri is an example of the former (obscuration of low clouds) 
while the area behind the dry line in north central Oklahoma illustrates 
the latter (refer to surface analysis. Figure 2-1). 

Results of previous studies have indicated that cumulus motion 
represents mean winds in the subcloud layer to within 5 degrees and 
1 ms~^ (Fujita et al., 1975 and Hasler et al., 1975). More recently, 
Suchman and Martin (1976) explored the accuracy and representativeness 
of tracer winds in the GATE* area. They found that ship winds used as 
"ground" truth differed from satellite winds by less than 3 ms Using 
30-minute interval imagery (as compared to 5 minute in this study) , they 


*£lobal Atmospheric Research Program Atlantic Tropical Experiment 


Table 3-1 




, 1-3 12 34 10.95 16 50 

3-4 14.80 9.77 17 73 


2-4 13 36 9.53 16 41 

1-4 13 16 10.56 16 87 


Location and Motion 


Direction 

(degrees] 


Latitude 

(“N) 


Longitude 

(°W) 


225 

33 

IMAGEl 

232 

08 

IMAGE2 

228 

42 

IMAGE3 

236 

57 

IMAGE4 


35.0503 96.2058 
35 0841 96.1640 
35 1090 96.1250 
35 1351 96.0767 


234 50 


CD 


231 26 



Fxgure 3-3 Satellxte derived wind fxeld at 1800 GMT. 
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Figure 3-4 Satellite derived wind field at 2100 GMI. 
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found a reproducibility of 1.3 ms ^ for cumulus level winds. Cumulus 
tracked winds would seem to better indicate boundary layer motion than 
do corresponding surface winds because: 

1. a 20-minute "average" of mean flow is obtained, 

2 surface winds are influenced by topography and surface stress. 
Doswell C1976) recognized the importance of using "filtered" surface 
data in computations of low-level moisture divergence. 

The wind fields of Figures 3-3 and 3-4 were then objectively 
analyzed by the cubic spline technique. This procedure was initially 
used by Fritsch (1971) in analyses of two-dimensional fields The 
obj* actively analyzed fields at 1800 and 2100 GMT are shown in 
Figures 3-5 and 3-6 respectively, with the vectors scaled as in 
Figure 3-3 For the 1800 GMT field, the grid interval is 0.5 degrees, 
while for the 2100 GMT 'field it is 0.4 degrees. This was done to mini- 
mize computer calculations in the 1800 field, due to its larger area 
and higher density of wind vectors For the 2100 wind field, the 
spatial distribution of the data indicate that the validity of the 
gridded winds would be most questionable over the northwestern and 
southeastern comers of the domain; however, these regions were not 
affected by storms. 

3.3 Moisture Convergence 

Several researchers, including Hudson (1971) and Sasaki (1973) 

computed surface moisture convergence and found good agreement between 

those patterns and reports of severe weather several hours later (see 

-3 -1 -1 

Figure 3-7a and b) Intense moisture convergence (~10 g kg s ) is 
one condition found favorable for tornadic storm development (Sasaki, 
1975) . Previous work by this author and others (Negri et al , 1971 
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Figure 3-5 Objectively analyzed wind field at 1800 GMT. 
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Figure 3-6 Objectively analyzed wind field at 2100 GMT. 
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a) Surface moisture convergence 

(10“^g kg“^s“^) 19 April 1972 
1800 CST. 


b) NWS radar summary 19 April 1972 
1945 GST. 


Figure 3-7 (from Sasaki, 1973). 
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and Maddox et al., 1977) showed that identifiable maxima in the moisture 
convergence field may precede development of intense convection by as 
much as 2-3 hours. 

3.3.1 Moisture convergence using satellite derived precipitable water 
On 24 April 1975, Vertical Temperature Profile Radiometer CVTPR) 
data was taken by the polar orbiting NOAA-4 sat'ellxte. This uses the 
15 ym CO 2 channels for temperature retrieval and the rotational water 
vapor absorption band to infer precipitable water (PW) values (Hillger 
and Vender Haar, 1976) This section explores the possibility of 
combining the satellite derived wind and moisture fields to simulate 
data that will become available with the launch of the VISSR* Atmo- 
spheric Sounder (VAS) on a geosynchronous satellite in the 1980 's. A 
mesoscaie moisture analysis is combined with the satellite derived winds 
to allow for the computation of moisture convergence in the boundary 
layer for the pre-storm environment. 

The moisture analysis at 1616 GMT was available in terms of total 
precipitable water values at a resolution of -70 km The conventional 
radiosonde network would provide similar data at a resolution of about 
400 km Figure 3-8 shows the derived PW field. An iterative retrieval 

r 

program (Hillger and Vender Haar, 1977) retrieves temperatures (and 
moisture) in mostly clear columns of the atmosphere High clouds 

t 

(above 500 mb) in central Kansas caused the PW values to be too low. 
Otherwise, the individual values should be accurate to within ±0.5 cm 
of precipitable water 

This moisture field was then objectively analyzed in the manner of 
the 1800 GMT winds and the horizontal moisture divergence (MD) was 

*Visible and Infrared ^can Radiometer 




Figure 3-8 Satellite derived precipitable water field (cm) 

at 1615 GOT (from Hillger and Vender Haar, 1977). 
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calculated at each grid point. 

MD = V • (PW V) (3-1) 

where PW is the precipitable water (cm) and V- is the velocity 
vector The significance of this term can be best understood by using 
the vector identity: 

V • (MV) = V ♦ VM + M(7*V) (3-2) 

where M is a moisture dependent variable (mixing ratio, precipitable 
water, dew point) and V is the velocity. 

The expression tends to have large negative values under conditions 
of 

1. advection of moist air [V • VM < 0] and 

2 horizontal convergence [M(V*V) < 0] . 

Moisture convergence is thus a measure of both the inflow of moist air 
and of the low-level convergence needed to create upward vertical motion 
The moisture divergence field is illustrated in Figure 3-9 Only nega- 
tive values (moisture convergence) have been contoured. The maximum 

-5 -1 

moisture convergence was ~30 x 10 cm H 2 O s centered in eastern 
Kansas. Local maxima are also found in southern Missouri, western 
Oklahoma, central Texas, and northern Louisiana. The radar summary at 
1835 GMT (Figure 3-10) shows some convective activity occurring in 
Kansas, however, these storms were not severe 

There is a rather large time discrepancy between the moisture field 
(1600 GMT), the wind field (1800) and the time of initial storm genesis 
(2200) Because of this, emphasis in this section has been placed on 
the description of a technique, rather than on demonstrating the 
possible correlation between moisture convergence and the occurrence of 
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severe local storms The next section will use the technique to couple 
the 2100 GMT derived winds with corresponding surface mixing ratios to 
determine possible correlations. It has been stressed here, however, 
that using data from satellite sensors alone, wind and moisture fields 
can be constructed and combined to give mesoscale horizontal moisture 
fltix information. A proposed vertical temperature sounder, placed on a 
geosynchronous satellite, will provide a vital tool for the implementa- 
tion and evaluation of this technique, notably as a real time "now-cast" 
forecasting aid. 

3.3.2 Moisture convergence from surface mixing ratios 

Surface mixing ratios were computed from the 2100 GMT surface 
observations using the definitions of mixing ratio and Teten's empirical 
formula for vapor pressure: 


q 


s 


0.622 


e^(T) 

p - 


(3-3) 


e 

s 


(T) - 610.78 exp[' 


17. 269 (T - 273.16). 
X - 35.86 ^ 


where: 


q is the saturation mixing ratio 
s 

is the saturation vapor pressure 
P IS the (surface) pressure 


(3-4) 


The saturation vapor pressure (e ) is in Pascals (Pa) when the tempera- 

s 

ture IS in ®K. The mixing ratio (q) is then computed using the dew 
point temperature (T^) . The resultant field of V • (qV) is presented 
in Figure 3-11. Units are g kg ^ x 10 ^ Comparison with Figure 3-12 
shows that during the following 2 hours storms developed in most regions 
characterized by large values of moisture convergence at 2100 GMT. The 
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Figure 3-11 


2100 GI^ moisture divergence field using 
satellite derived winds and surface mixing 
ratios. Only negative values Cii^oisture 
convergence) have been contoured Units 


are g kg“^s~^ x lO"^ 
50 X 10“^ g kg“^s“^. 


and contour interval is 
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Figure 3-12 


Severe storm locatxons at 2300 GMT. 
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mam exception was the lack o£ activity within the convergence region in 

central Kansas A reexamination of the original wind field (Figure 3-4) 

shows that almost no wind vectors are foxmd in that area, so the 

validity of the gridded winds there are questionable Maximum values 

-3 -1 -1 

of moisture convergence (-2 2 x 10 g kg s ) are similar to values 
found by Sasaki (1975) to be associated with severe thunderstorm occur- 
rences. The Neosho, Missouri tomadic storm began to develop within an 
hour of the satellite measurements and within the area of highest 
moisture convergence The same field using only 2100 GMT surface winds 
is shown in Figure 3-13. The overall pattern of the two fields is 
similar, however, the surface based analysis indicates moisture 
divergence over much of south and southcentral Oklahoma This is an 
important difference since several severe storms developed across that 
region There are also areas of moisture convergence where no severe 
weather occurred An improved correspondence between maximum moisture 
convergence and subsequent severe storm genesis was characterized by 
the satellite field Elimination of local effects on surface winds and 
a stronger flow at cloud base are postulated as the cause 

3.4 Dynamic Parameters 
3.4 1 The divergence term 

Wilson (1976) has computed large-scale vertical motion fields using 
the Atmospheric Variability Experiment (AVE IV) sounding data collected 
on 24 April. His computations showed that the large-scale vertical 
motions did not play a major role in organizing the convection on this 
day since the storms developed in regions of diagnosed synoptic scale 
subsidence or of only slight upward motion (refer to Figure 2-6) There- 
fore, these storms were induced by mesoscale processes. Endlich and 



34 



Figure 3-13 2100 GMT moisture divergence field using 

observed surface winds and mixing ratios 
Units and contour interval are the same 
as in Figure 3-11 
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Mancuso (1968) showed that the term for vorticity production in the 
vorticity equation was a good objective indicator of severe storm 
activity The vorticity equation in Cartesian coordinates is; 


A(c+£) - -(?+£) (|J 


.8w ^ ^ 9U x , 1 ^3p ^ _ 9£. 3P \ 

^3x 3z 9y 32 *^ 


(3-5) 


where ? is the relative (parcel) vorticity 
f IS the Coriolis (earth) vorticity 
u,v,w are the components of the velocity vector 
p is the density 

Equation (3-5) states that the rate of change of the absolute 
vorticity following the motion is given by the sum of the three terms 
on the right, called the divergence term, the tilting term and the 
solenoidal term respectively. Using Doppler wind information, Ray et al 
(1976) found the tilting term to be an order of magnitude smaller than 
the divergence team. Scale analysis (Holton, 1972) shows the solenoidal 
term can also be neglected 

From the definition of relative vorticity (i;) and divergence (D) 
in finite difference form 



f ^ f) 


(3-8) 
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For mesoscale systems ? cannot be neglected with respect to f as is 

often done in synoptic scale computations. The field of ’’vorticity 

-2 -9 

production" is shown in Figure 3-14. Units are s x 10 and only 

positive values are contoured Comparison with Figure 3-12 shows a 

significant correspondence exists between areas of high mesoscale 

vorticity production and the subsequent occurrence of severe local 

-9 -2 

storms It appears that values on the order of 15 x 10 s are 
necessary and sufficient to force the occurrence of severe weather 
3.4 2 Relative vorticity 

The field of relative vorticity (5) alone can provide some clue as 
to severe storm genesis, as well as provide an insight into the sensi- 
tivity of the data set to the density and number of wind vectors 
Figure 3-15 is the 2100 GMT relative vorticity computed using finite 
differencing on the satellite derived wind field The vorticity is 

V 

concentrated in a narrow band along the front and dry line as would be 

-4 -1 

expected. Maximum values of 2 x 10 s are of the order of the earth’s 
vorticity at that latitude. A subjective modification was made to the 
original data set such that eleven vectors near the rear of the fronts 
were arbitrarily deleted The modified field is shown in Figure 3-16. 
The results differ in both general pattern and by as much as two orders 
of magnitude in some regions. The modification was made to maximize 
differences, however, the ability to track small cumulus is dependent 
on real weather systems. It is often impossible to find cumulus clouds 
behind a dry line, and low-level clouds are often obscured by layers of 
middle and high clouds to the north of a surface front. The Neosho, 
Missouri tomadic storm developed in an area of a local maximum of 
relative vorticity. While the concentration of vorticity for tornado 
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Fxgure 3-14 "Vortxcity productxon" at 2100 GMT. 

Only positive values have been 

-2 -9 

contoured. Unxts are -s x 10 and 

-9 -2 

contour xnterval xs 5 x 10 s 
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Figure 3-15 Relatxve vortxcxty computed from satellxte 
derived wxuds at 2100 GMT» Only posxtive 
values have been contoured. Units are 
-1 -4 

s X 10 and contour interval xs 
5 X 10"^ s“^. 
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27tt-01 



Figure 3-16 Relative vorticity computed from the 
modified satellite derived winds. 
Details are as in Figure 3-15. 
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genesis is on a scale smaller than the mesoscale, the satellite derived 
field indicates a high vorticity environment was present at the outset. 

3.5 Summary of Wind and Moisture Fields 

This chapter has detailed some of the computations that can be made 
using satellite derived wind fields. Using data from satellite sensors 
alone, horizontal moisture flux information can be obtained Because 
of the early NOAA-4 pass time, and a 2-hour discontinuity between the 
wind and moisture fields, no inferences were made to the correlation of 
low-level moisture convergence and the occurrence of severe local storms 
at 1800 GMT Using 2100 GMT surface mixing ratios as the moisture 
parameter, however, subsequent severe storm development occurred within 
several of the regions of diagnosed high moisture convergence. Satel- 
lite moisture convergence patterns specified convective patterns that 
developed better than did the same field computed using only surface 
observations. This suggests that satellite data may be useful in 
operational severe storm forecasting. Several dynamic parameters 
(vorticity production, relative vorticity) can also be used as good 
indicators of severe storm potential Table 3-2 summarizes the results 
of this and previous studies. 

3.5,1 Limitations 

The author is well aware of the limitations of the technique 
described. One such limitation is the fact that the satellite winds 
are not exactly located in the vertical. Schaefer (1973) pointed out 
that when divergence is computed on a non-horizontal surface, additional 
terms appear in the defining equation. Consider (x,y,z) as the spatial 
coordinates in a constant height system and (x,y,p) as the spatial 
coordinates in a constant pressure system. The formula for transferring 



Table 3-2. Previous Studies of Severe Weather -Related Parameters 


Investigator 

Scale 

Divergence 

(s'^) 

Relative 

Vorticity 

(s“^) 

Moisture Divergence 

Vorticity 

Production 

(s ) 

Endlich and 
Mancuso , 1968 

synoptic 



-20 X 10"^[°C s"^] 

30 X 10“^*^ 

Viezee et al., 
1972 

synoptic 


2 X 10“^ 



Sasaki, 1973 

mesoscale 

-5 X 10"^ 


-40 X 10"^[gkg"^s‘^] 


Hougton and 
Wilson, 1975 

mesoscale 

-5 X 10“^ 




Ostby, 1975 

mesoscale 



-150 X 10"^[gkg"^s“^] 


Suchman and 
Martin, 1976 

mesoscale 

-6 X 10"^ 




Peslen, 1977 

mesoscale 

-8 X lO"^ 

"“5 "A 

10 ^ - 10 



Maddox et al., 
1977* 

mesoscale 


-4 

2 X 10 

-200 X 10“^[gkg“^s“^] 


Negri et al. , 
1977* 

mesoscale 

-10 X 10“^ 


-20 X 10“^ [cm H^O s"^] 


Present study 

mesoscale 

-10 X 10“^ 

2 X lO"^ 

-200 X 10"^[gkg~V^] 

-10 

200 X 10 


*Same data set as present study 
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the divergence from one system to another is given by: 

D = = (7j,.V)p - § . (7„z)p C3-9) 

This states that the divergence on a surface of constant height equals 

the difference between the divergence on a surface of constant pressure 

and the product of the vertical wind shear and the slope of the height 

surfaces For this case, it was assumed that the satellite winds were 

on a horizontal surface, even though rawinsonde data indicated that the 

LCL (lifted condensation level, an estimate of cloud base) ranged from 

910 to 800 mb across the region in question. For small vertical wind 

-5 -1 

shears and gentle slopes (last term in (3-9) « 10 s ) the simplifying 
assumptions remain valid. Other limitations include the fact that the 
1800 GMT wind fields were multiplied by an estimation of the total 
atmospheric water content, not just boundary layer moisture. Total 
precipitable water content is likewise not retrieved perfectly It was 
also shown that the derived fields are sensitive to the density and 
spatial configuration of the original wind fields, and most likely 
sensitive as well to the objective analysis scheme used. Satellite 
derived winds lend themselves to processing by multi-variate analysis, 
which would simultaneously analyze satellite winds, surface and upper 
air data, and other data in the best possible way to reduce error. 



4.0 SATELLITE MONITORING OF SEVERE LOCAL STOEMS 


4.1 Previous Studies 

It was shown in the last section how monitoring of the pre-storm 
environment through the use of satellite derived convergence patterns 
could give an indication of subsequent severe storm genesis areas. 
Equally, or perhaps more important, is the continued observation of 
these areas by a geosynchronous satellite having high temporal frequency. 
Satellite measurements in the atmospheric "window" (10-12 pra) can dif- 
ferentiate between convective elements of differing duration and inten- 
sity. One may then deduce the intensity of this convection and then 
(hopefully) be able to discriminate between severe and non-severe 
elements . 

Previous attempts using satellite infrared (IR) data have largely 
been qualitative and have relied on the use of enhanced IR imagery. 

Purdom (1976a) used rapid interval high resolution GOES imagery to view 
tornadic storms from both an earth relative frame and a thunderstorm 
relative frame. He concluded that a relative low-level south to south- 
east cloud flow into a severe storm is important in that stom's ability 
to produce a tornado. Purdom (1976b) has also observed the presence of 
a variety of trigger mechanisms and interactions that range in scale 
from the synoptic to the mesoscale. Among such mechanisms that exert a 
strong influence on the initiation and maintenance of thunderstorms are 
convective lines, convective cloud mergers, and convective line inter- 
sections, all visible on satellite imagery. Oliver and Scofield (1976) 
and Scofield and Oliver (1977) present schemes for estimating convective 
rainfall from satellite visible and IR imagery. The major assumptions 
of their scheme are summarized in Table 4-1. While this present study 
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Table 4-1. Convective Raxnfall Estimation Scheme 
(from Scofield and Oliver, 1977) 

1. Bright clouds in the visible imagery produce more rainfall 
than darker clouds. 

2. Bright clouds in the visible and clouds with cold tops in the 
IR imagery which are expanding in areal coverage (in early and mature 
stages of development) produce more rainfall than those not expanding. 

3. Decaying clouds produce little or no rainfall. 

4. Clouds with cold tops in the IR imagery produce more rainfall 
than those with warmer tops. 

5. Clouds with cold tops that are becoming warmer produce little 
or no rainfall. 

6. Merging o£ cumulonimbus (Cb) clouds increases the rainfall 
rate of the merging clouds. 

7. Most of the significant rainfall occurs in the upwind (at 
anvil level) portion of a convective system. The highest and coldest 
clouds form where the thunderstorms are most vigorous and the rain 
heaviest. These cold clouds get thinner downwind and look warmer in IR 
imagery as the anvil material blows away from its origin. 
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does not address itself specifically to the determination of convective 
rainfall amount, several of the assumptions of Table 4-1 can be applied 
to the detection of severe thunderstorms. 

Other authors took the next step toward a more quantitative approach 
to the remote sensing of convection. Am (1975) used the Optical Data 
Digitizer and Display System (OD ) to digitize film loops from the third 
Applications Technology Satellite (ATS**3) and from the second synchro- 
nous meteorological satellite (SMS-2) He analyzed visible channel 
radiances ("brightness”) from the anvils of severe thunderstorms in an 
attempt to discern signatures which would describe storm severity. 

Among the parameters found useful were anvil area divergence (1/A dA/dt) , 

2 2 

area speed gradients (d A/dt ), and volume flux (dVo/dt). Griffith 
et al. (1976) used a densitometer to examine ATS visible images. They 
found that daily rainfall in the tropics could be estimated from cloud 
brightness, a technique which is currently being tested in the 
mid-latitudes . 

Perhaps the most quantitative method of observing severe storm 
characteristics is through the use of the original satellite digital 
data itself. Auvine and Anderson (1972) obtained upper tropospheric 
fluxes and divergence for several storms from measurements of anvil 
area expansion. They presented two models based on the point source 
concept of a growing cumulonimbus anvil. These models are reproduced 
in Figure 4-1, and conform well XTith actual observations of anvil shape 
and expansion. Adler and Fenn (1976, 1977) have developed techniques 
to use high frequency SMS infrared data to help detect and monitor 
severe thunderstorms. This study will use some of their techniques in 
an attempt to discriminate between severe and non-severe storms from the 
vantage point of a satellite platform. Adler and Fenn showed that three 




a) A model of airflow wxthiti the thunderstorm 
anvil assuming calm environmental winds. 

/ 


b) The material outline of a cumulonimbus 
anvil assuming its growth to be the 
result of a point source in a uniform 
stream . 


Figure 4-1 Thunderstorm anvil models (from Auvine and Anderson, 1972). 
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simple characteristics seemed to be correlated with the occurrence of 
severe weather. These are duration, minimum cloud top temperature, and 
cold area rate of increase. The following sections will examine these 
characteristics for the storms of 24 April 1975, after a brief look at 
the theory of remote sensing. 

4.2 Radiative Transfer Equation 

In the spectral interval 10-12 pm, the clear atmosphere is 
essentially a non-scattering, non-absorbing medium. Thus, the equation 
of radiative transfer reduces to: 


N(X) = B(X, T)xe 


(4-1) 


where N(X) is the spectral radiance observed at the top of the 
atmosphere in a spectral interval dX around mean 
wavelength X 

T IS the transmittance at wavelength X of any optically 
absorbing gas (assumed t = 1 for X = 10-12 ym) 
e IS the surface or cloud top emissivity (assumed e = 1) . 
The Planck blackbody radiance B(X, T) at mean wavelength X and 
temperature T is 


B(X, T) 


C-/XT 

e - 1 


(4-2) 


Thus, the temperature at which the top of a cloud radiates may be 
expressed as an effective blackbody temperature (Tg^) which can be con- 
tinuously monitored by an infrared radiometer on a geosynchronous 


satellite. 
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4.3 Data and Methods 

SMS digital visible and infrared data were obtained for the period 
1800-0200 GMT on 24-25 April 1975. During this period the satellite 
was in a limited scan mode, wherein a small north-south segment is 
observed on a frequent basis (nominally five minutes) Data coverage 
on this time scale is necessary for severe storm monitoring when one 
considers the lifetime of most convective elements is less than 30 
minutes. The clouds examined in this study are shown in the 2242 GMT 
infrared image of Figure 4-2. Clouds designated by names or numbers 
became severe storms, while the non-severe storm clouds are lettered. 
Visible channel data were not quantitatively examined due to low solar 
zenith angles during storm lifetimes. 

The spatial resolution of the IR instrument is approximately 32 km‘ 
at nadir. Due to sampling constraints, however, IR "pixels" or scan 
spots were sampled twice in the east-west direction. Thus, each IR 

2 

element referred to in later sections represents a mean area of ~90 km 
at 40® north latitude. Infrared growth rates and temperatures are 
examined in an attempt to. 

1. discriminate between severe and non-severe elements, and 

2* link signatures in the IR data to the occurrence of severe 
weather on the ground • 

This study will follow techniques initially used by Adler and Fenn 
(1976, 1977) to portray SMS infrared data in a temporal sequence. This 
was done to promote consistency in examining digital data as well as to 
facilitate the direct comparison of results. Starting at the cold end 
of the temperature spectrum, a cumulative histogram of data at or below 
a given temperature is compiled. That is, the number of points 
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ME02 004413.11 IRl T/L: 0170/ 1^00 24 hPR7S 2242Z 0000 


CLOUD 2 


NEOSHO^^ 
CLOUD L-«“ 


CLOUD A 


CLOUD B 


CLOUD C 


WEWOKA-^ 


CL^UD 


CLOUD C 


Figure 4-2 SMS-2 infrared image at 2242 GMT. 

Cloiad examined in this study are labelled. 
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with blackbody temperature T„ < T. is obtained. The temperatures 
(T^) range from the coldest temperature present to the warmest tempera- 
ture, T^, which is unambiguously part of the cloud. This temperature 
(T^) is defined such that its contour is far enough within the cloud to 
avoid the effects of lower emissivity at the anvil edge. (An exception 
to this procedure for weak convective clouds will be investigated later.) 
When the warmer temperature contours near the edge of the anvil merge 
with the anvils of other storms, the warmest valid temperature (T^) 
decreases. Late in the lifetimes of the storms, only the very cold 
areas are followed and thus only a small range of T exi„st. Due to 
the method of data reduction utilized, the T.j, values represent a 
mean value of temperature over a 4“C interval. When applied in a 
consistent manner, this will not introduce significant error into the 
derived rates of growth. Rather than display the contour of every T„_, 
representative temperatures were chosen to portray a wide range of anvil 
area. Particular attention, however, was given to the very cold 
portions of the cloud top (206-214 K), as they are often indicators of 
intense convection below. 

4.4 Results for Severe Storms 

Four of the clouds examined on this day produced tornadoes and/or 
large hailstones. Three of those storms will now be examined in detail: 
the Wewoka hailstorm, the Neosho tomadic storm and cloud no. 1, a hail 
producer (refer to Figure 4-2). 

4.4.1 Wewoka 

Figure 4-3 is the growth rate diagram for the cloud that produced 
baseball-sized hail over the town of Wewoka, Oklahoma. The graph dis- 
plays the relationship between the number of data points, N^, and time. 


NUMBER (N) 



Figure 4-3 Thunderstorm growth rate diagram for the Wewoka hailstorm. N is the number 
of data points in the element with blackbody temperature _< T^, where 

IS as identified in the diagram. 
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Again, is the number of IR pixels in the cloud top with blackbody 
temperature T „ less than or equal to T . For example, in 

DD 1 

2 

Figure 4-3, the 222K isotherm, which contained 7 pixels (-630 km) at 

2 

2240 GMT, expanded to 21 pixels (-1890 km ) in just five minutes. Be- 
cause the ordinate of the graph is logarithmic, the slopes of the lines 
are proportional to ^ dN/dt Adler and Fenn (1977) found that in seven 
of eight cases studied, the first report of severe weather occurred 
during, or just after, a rapid expansion of cold areas. It is worth 
noting that the Wewoka storm's coldest temperature (-67°C to -63®C) did 
not appear until two hours after the occurrence of hail Thus, while 
the temperature at cloud top may indicate its convective potential in 
terms of rainfall (Scofield and Oliver, 1977), it is not necessarily a 
precursor to storm severity, i.e., tornadoes and large hail. Figure 4-4 
illustrates the temporal changes in divergence and rate of decrease of 
cloud top minimum temperature for the Weiroka storm. The bottom portion 
represents the slopes of the family of curves of Figure 4-3. For 
example, the change in divergence of the 222K isotherm is found by 

following the open triangular symbols. Divergence values of 
-3 -1 

1-5 X 10 s are found in the pre-hail period, while much slower 
rates of 5-10 x 10 ^ s ^ are observed after the occurrence of hail. 

This order of magnitude difference is physically related to the vertical 
motion occurring under the cirrus anvil. The rate of change of minimum 
blackbody temperature (Figure 4-4, top) provides little insight as to 
the time of hail occurrence. An initial period of rapid decrease is 
followed by a period where temperature changes little with tune. Mini- 
mum blackbody temperature was thus a poor indicator of the time of severe 
weather in this case; however, it will be shown that it is a good 
indicator of severe versus non-severe elements 
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Fxgure 4-4 Temporal change xn dxvergence and rate of decrease of cloud top mxnxmum 
temperature (top) for the Wewoka storm. 
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4.4.2 Neosho 

As viewed from the satellite, the Neosho storm was much less 
spectacular that the neighboring Wewoka storm. The growth of the Neosho 
storm, considerably less than, the explosive growth of Wewoka, is shown 
in Figure 4-5. In this case, the rate of growth did not specify storm 
severity. The storm spawned tornadoes at Miami, Oklahoma, at 0010 GMT 
and at Neosho, Missouri, at 0040, nearly two hours after its inception. 
Studies by Fujita (1973) indicate a preference for the formation of 
stronger tornadoes later in the day. The analysis is also hampered by 
a lack of data between 2330 and 0000 GMT. While both reports of torna- 
does occur during or following rapid cold area increases, the observa- 
tions rest weakly on the basis of only several data points. The 
variation in divergence, and in rate of change of minimum temperature 
(Figure 4-6) also show very little in the way of identifiable signatures. 
Active convection near the time of the tornadoes was most likely obscured 
by the extensive cirrus anvil. Overshooting tops were apparent in the 
visible imagery, however. Fujita (1973) proposed a mechanism of tornado 
formation whereby the collapse of the overshooting top, due to precipita- 
tion overloading, induces a twisting doxmdraft inside of a rotating 
storm. Overshooting tops are not apparent in the IR channel due to its 
coarser resolution and finite response time. An overshooting tower 
rapidly becomes much colder that its environment as it overshoots and 
penetrates the tropopause. However, the storm tops entrain environmental 
air as they "boil" above the anvil. Because the environmental air tem- 
perature changes little with height, the net effect is the smoothing or 
masking of the presence of the very cold temperature within the cloud 


top. 
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Figure 4-6 shows the divergence and temperature change as a function 

of time for the Neosho storm. Near the times of the tornadoes, cloud 

-3 -1 

top divergences were small, (0.1-1. Ox 10 s ) and fairly large 
scatter existed. The storm top had apparently reached its minimum 
temperature also. The above points to the need for a bispectral or even 
multispectral method of observing severe storms. A microwave radiometer, 
with proper polarization, would be able to remotely sense the precipita- 
tion size particles below the anvil (Weinman and Guetter, 1977) . Com- 
bined with visible channel radiances to detect overshoots and IR channel 
radiances for temperature, a nearly complete description of a storm from 
satellite sensors could be inferred 
4 4.3 Cloud no. 1 

The growth rates of this hail producing storm are displayed in 

Figure 4-7. Though initially not as steep as those of the Wewoka storm, 

-3 -1 

normalized rates of growth remain greater than 3 x 10 s for several 

hours after the inception of the cloud. First reports of hail at 2325 

-3 -1 

GMT occurred during phases of growth as large as 5 x 10 s . Of 
particular interest are the merging of three isotherms at 2325 with the 
anvil of the Wewoka storm. Scofield and Oliver (1977) and Purdom (1976b) 
discuss the importance of convective element mergers to the enhancement 
of convective rainfall and thunderstorm maintenance respectively. 

4.5 Comparison with Non-Severe Storms 

Most of the precipitating clouds examined in this study became 
severe at some point in their lifetime. There were, however, several 
clouds that formed at or near the times of those that did become severe. 
In Figure 4-2, clouds lettered A through E were all non-severe, with 
severity determined by reports published in Storm Data . Growth rates 
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for three non-severe storms are presented in Figure 4-8. Differences 
between this figure and those for the severe storms are fairly obvious. 
First, the minimum blackbody temperature (Tgg) used to define these 
clouds had to be increased from 244 K (-29“C) to 256 K (-IT^C) This 
IS probably an vmderestiraate of the true temperature due to non -unity 
emissivities at cloud edges with correspondingly "warmer" Tgg* A. 
further problem might be one found by Negri et al (1976) They found 
that the IR sensor did not correctly respond to clouds of small horizon- 
tal extent The satellite has a rotation period of 0.6 s and only views 
the earth for five percent of that time. One scan is composed of 3822 

IR pixels and thus the sensor will only "see" a cloud of 30 km width 

(7 pixels) for 5 x 10 ^ s. This is approximately the response time 
(55 ysec) of the IR sensor and thus the correct response to the cold 
cloud radiance is not realized for clouds of small horizontal extent. 

Several other apparent differences exist between severe and 
non-severe growth rates. In Figure 4-8, the duration of the clouds, 

as well as the number of pixels defining the clouds, are much less than 

-3 -1 

those of the severe storms. Growth rates rarely exceed 1 x 10 s , 
and in fact one element (Cloud D) had a negative rate of growth. Real 

t 

txme monitorxng of growth rates on a frequent (5^15 mxnute) basis could 

xdentxfy clouds of probable severe potentxal early xn thexr Ixfetxme. 

Monxtorxng of mxnxmum blackbody temperature as well as growth 

I>D 

rates can also provxde an xnsxght xnto storm severxty. Fxgure 4-9 xs a 

representatxon of mxnxmum T wxth txme for two severe and two non- 

severe storms. The severe storms exhibited a much greater decrease in 

minimum than did the neighboring clouds. Minimum temperatures 

BB 

reached were 10-15°C colder for the severe storms. As was mentioned 
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Fxgure 4-9 Temporal variation of minimum blackbody temperature (T ) for 
two severe and two non~severe storms. 



62 


previously, the time of occurrence of the severe event cannot be 
inferred from these temperature patterns 

4.6 Results for All Storms 

On the basis of the above parameters a graphical summary of all 
clouds can be constructed which relates the severity of a storm to 
several key parameters. Adler and Fenn (1977) have done this for 
several storms that took place on 6 May 197S. Their results are repro- 
duced in Figure 4-10 The graph shows the relation among cold area 
rates of growth, cloud top relative temperature, and the occurrence of 
severe weather The quantity is the cloud top relative temp- 

erature, where T^, is the minimum temperature achieved during the 
lifetime of a cloud element, and is the minimum temperature of 

the coldest element during the entire lifetime of the cloud system 
For a cloud that grows as cold or colder than the other clouds in the 
area, (T^- T,,,..J will be zero. The values of ^ dN/dt represent the 
maximum values, irrespective of temperature, for the very early stages 

of cloud growth. Their results show good agreement between nomialized 

-3 -1 

rates of growth greater than 3 x 10 s and storm severity, and also 
between low values of cloud top relative temperature and storm severity 
This type of representation is fine for an after-the-fact analysis, but 
It cannot be used in real-time decision making because’ the minimum 
temperatures T^ and not known until after the occurrence 

(or non-occurrence) of severe weather Thus, the technique was modified 
in this study for use as a real-time forecasting aid, only data for the 
first half hour of cloud lifetimes were used This would provide an 
operationally significant time interval for decision making and 
dissemination of information Results are presented in Figure 4-11. 
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Figure 4-10 Relation of cold area rate of growth to cloud top relative temperature 
and the occurrence of severe weather (from Adler and Fenn, 1977). 
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Fxgure 4-11 Results from this study using only the first 30 minutes of cloud lifetime 
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The quantities CTq-Tj^jjj) and ^ dN/dt are still defined in the manner 
of Adler and Fenn, but are indicative of the very early stages of con- 
vective development. All non-severe storms had initial growth rates 

-3 -1 

less than 1.3 x 10 s , with correspondingly warmer minimum tempera- 
tures than the severe storms. For severe storms, normalized growth 
-3 -1 

rates of 2 0 x 10 s or greater were found It should be noted that 
sample size is very limited, and the need for a larger statistical 
sample from more case studies is needed. Despite' this, it appears as 
though the early detection of severe weather producing clouds can be 
determined by frequent infrared monitoring by a geostationary satellite. 

4.7 Summary of Infrared Growth Rates 

It appears from this study that rapid- interval SMS infrared digital 
data can be useful in the monitoring of severe local storms. Data of 
this temporal resolution (~5 minutes) are necessary for a "now-casting" 
capability for severe weather events. Several parameters may be cal- 
culated which are related to thunderstorm intensity and the occurrence 
of severe weather on the ground. In the early phase of cloud lifetime 
(first 30 minutes) criteria for a severe-non-severe decision should 
include: 

1. the normalized rate of growth (or cloud-top divergence), 

1 —3 

“ dN/dt, Maxxmum values greater than or equal to 2 x lO" were indxca- 

txve of subsequent severe weather* 

2 . relative cloud-top temperature . Cloud elements 

10-12 °C colder than nexghborxng clouds are suspect as potential severe 
storms, 

3. rate of change of mxnimum temperature. Cloud elements wxth 
decreases of 1.5^C/mxn xn their early stages were found to have become 


severe 
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4 duration and size. Non-severe cloud elements were found to be 

of smaller area and of shorter duration than severe elements. Larger 

and longer •=• lived clouds were of greater severe potential. 

Predicting the time of initial hail or tornado occurrence was realized 

with a lesser degree of success. Reports of large hail occurred after 

-3 “1 

several consecutive periods of rapid growth (1 - 5 x 10 s ) However, 
the most destructive storm, the Neosho tornado, provided no discemable 
signatures in growth rate or temperature pattern. This is attributable 
to the late occurrence of the tornado below the cirrus anvil and to the 
inadequate infrared resolution of overshooting tops Prior analysis by 
wind tracking, however, had indicated the high vorticity environment in 
which the Neosho storm developed. Results were in general agreement 
with those of Adler and Fenn, but were hindered by a lack of data prior 
to the tornado and by the small number of clouds used in the data base. 
No distinction was possible concerning the individual effects of wind 
shear and vertical velocity on the magnitude of the normalized rates of 
growth. It was assumed that the expansion of the cloud-top (due to 
generation of ice crystals) was due primarily to the intensity of the 
vertical motion. 

Higher quality data will hopefully be available in the 1980 's with 
the launch of StormSat, a 3-axis stabilized sounding and imaging system 
dedicated to meso scale phenomena. Beyond that, the Synchronous Earth 
Observatory Satellite (SEOS) will provide improved visible and infrared 
sounding combined with microwave imaging and sounding 



5,0 RADAR^-SATELLITE COMPARISON 


5.1 Areal Comparisons 

Radar data for the Wewoka hailstorm were available in the form of 
intensity contoured reflectivities (at zero degree elevation angle) 
from the WSR-57 radar at Oklahoma City (OKC) . No radar data were avail- 
able for the tornadic Neosho storm. The sequence of radar echoes 
between 2224 and 2329 GMT is shown in Figure 5-1. This interval repre- 
sents the inception and growth phase of the Wewoka storm and includes 
the time span when two-inch and larger-sized hail was reported on the 
ground (2300-2315 GMT). In Figure 5-1, the location of all echoes is 
shown with respect to the fixed location of the OKC radar (denoted by 
the dark x^s). The first echo developed south of the radar with sub- 
sequent echoes forming from southwest to northeast. The corresponding 
sequence of images as seen by the infrared (IR) sensor of the satellite 
IS illustrated in Figure 5-2. The truly explosive growth of this storm 
can be seen as the tremendous increase of IR area in slightly over an 
hour. The horizontal scale and position of true north are indicated on 
the figure, but differ from those of Figure 5-1. When scale factors are 
taken into account, the entire radar echo area corresponds in areal 
extent to only a small portion of the infrared cloud area. 

A more quantitative example of the radar echo growth of the Wewoka 

2 

storm IS portrayed in Figure 5-3, the radar echo area (km ) as a func- 
tion of time. The ordinate of the graph is logarithmic to facilitate a 
direct comparison with Figure 4-3, the infrared growth rate. Both the 
level 1 area (> 10 dbz) and the level 3 area (> 30 dbz) show initial 

periods of rapid growth followed by a slower growth after the occurrence 

-3 --1 

of hail. Initial growth rates are of order 5 x 10 s , virtually 
identical to those computed for the IR cloud area. The level 1 and 
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greater area was best related xn actual value and xn temporal varxatxon 
to the area encompassed by the 244K (-29“C) isotherm on the satellxte 
IR cloud. For a sample size of nxne data points, the correlatxon 
coeffxcxent of areas was 0.994. The level 3 and greater echo area, 
whxle well correlated with the temporal variatxon of the 210K (-63“C) 
isotherm, was of too small a magnxtude to correspond specxfxcally wxth 
any partxcular IR area. Thxs xs an xndxcatxon that xntense radar echoes 
of small horizontal extent cannot be detected as an IR sxgnature under- 
neath the expandxng cxrrus anvxl. Small, xntense echoes, as well as 

overshootxng domes, have diameters of 10-12 km or cross-sectional areas 

2 

of approximately 80 km . This is approximately the area of one IR data 
2 

point (90 km ) and thus the problem becomes one of resolution and 
sampling . 

5.2 Storm Top Comparisons 

The exact determination of cloud top height is of fundamental 
importance in weather forecasting and analysis . The vertical extent of 
a cloud IS a measure of the intensity of its upward vertical motion. 
Thus, indirectly, vertical cloud extent is related to low-level converg- 
ence, subsequent release of latent heat, and eventually to the intensity 
of rainfall and to the degree of storm severity. 

For the determination of cloud height from infrared radiances, we 
assume that a cloud emits radiation as a blackbody and that the cloud 
top temperature is the same as the environmental air temperature at that 
height. This was demonstrated by Pitts et al. (1975) from airborne 
radiometric measurements over severe storms. Particular problems can 
arise for small isolated clouds in a hot, dry environment due to radi- 
ometer response problems (Negri et al., 1976), but for large storms 
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extendxng to the tropopause the assumptions remain valxd. 

A direct comparison between radar and satellite heights was 
facilitated by the semi-xsolated nature of the storms involved. 

National Weather Service (NWS) radar summaries were compared to nearly 
simultaneous (maximum difference of 9 minutes) blackbody temperature 
fields in an attempt to quantitatively assess cloud heights. Results 
are presented in Figure 5-4, a comparison of radar indicated cloud 
heights with heights derived from minimum IR temperatures and the 
Monett environmental sounding. The measurements were divided into three 
arbitrary time intervals, to examine the effect of the phase of cloud 
development on the radar-satellite height comparisons. Below the tropo- 
pause (~42,000 ft) there is fair agreement between heights determined 
by the two methods of sensing, with radar tops consistently higher than 
those derived by satellite measurements. Bor storms growing above 
the tropopause, determination of cloud tops from IR dat'a alone was 
not possible. 

The apparent ambiguity between the two methods of cloud top 
determination can be resolved if several factors are considered. 

First, the "inexactness" of reported heights on the NWS radar sum- 
maries may account for some of the scatter in Figure 5-4. These 
heights represent instantaneous values of the echo top as seen by 
one or more radars during a 20-minute period prior to the actual 
report. Given the magnitude of vertical velocity in a severe storm, 
echo tops can change in the vertical several thousand feet in just a 
few minutes. Smith and Reynolds (1976) have also found root mean dif- 
ferences of 2.53 km in studies of IR/radar height pairs. 
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Infrared resolution, and response may also impose limitatxons in 

inferring cloud heights. Overshooting domes of 10-12 km diameter will 

2 

have cross sectional areas of approximately 80 km . The field of view 

2 

(FOV) of the IR sensor at 35®N is approximately 90 km so that an 
overshooting turret would have to be exactly centered in the FOV to 
return the correct radiometric response. Figure 5-5 shows a conceptual 
model illustrating this resolution problem. A cloud turret penetrating 
the tropopause will rapidly become much colder than its environment 
by cooling moist adiabatically. However, the small turrets are entraining 
environmental air whose temperature changes little with height because 
of the nearly isothermal lapse rate above the tropopause. The net ef- 
fect IS to smooth out the presence of the overshooting dome as viewed 
by the satellite IR sensor. 

The finite response time of the IR sensor must also be considered 

when viewing small, isolated turrets above the anvil. We can relate 

the observed temperature (T ) to the effective blackbody temperature 

ni0^s 

of the cloud by 


T == T 
meas eff 


- A dT/dN 


(5-1) 


where A is the lag coeffxcient (response) of the sensor and dT/dN 
IS the temperature gradient along the dxrection of the scan. 

A scan consists of 3822 pxxels (1911 FOV’s) completed whxle the 
sensor “sees” the earth for 5% of xts rotation perxod (0.6 s) . This 
allows each FOV to be vxewed for 1.57 x 10 ^ s . Since the response 
time of the sensor is 55 x 10 ^ s, 3-4 FOV’s are needed to realxze a 
correct response at the satellite. Overshootxng domes rarely provide 
this area xn the xnf rated. 




76 


6.0 CONCLUSIONS 

The approach taken in this study o£ severe local storms was 
essentially two-fold. The first technique utilized toward better pre- 
diction of severe storm genesis was the monitoring of the pre-storra 
environment through the use of satellite derived low-level wind fields. 

It was found that with the possible exception of not being exactly 
located in the vertical, the satellite wind fields provided a better 
representation of boundary layer flow then did corresponding surface 
winds Wind fields coupled with surface mixing ratios as the represen- 
tative moisture parameter allowed for the computation of moisture 

divergence in the pre- storm environment Areas with moisture conver- 

-1 -1 

gence values of -1.5 x 10 g kg s or greater were found to be sub- 
sequent genesis areas of severe storms. Similar calculations using 
observed surface winds did not specify convective patterns as accurately. 
An equally representative moisture variable, total precipitable water, 
can be derived from vertical sounding of the (clear) atmosphere and 
combined with satellite wind fields. Thus, when a geosynchronous 
atmospheric sounder becomes available, moisture divergence fields may 
be determined solely from the vantage point of a satellite platform 
As was demonstrated, simultaneous wind and moisture fields are desired 
if correlations between moisture convergence and severe storm genesis 
are to be considered valid. 

The dynamic parameters of boundary layer "vorticity production" 
and relative vorticity also provided insight into the developmental 
nature of severe storms. When a region had both low-level convergence 
and large, positive relative vorticity, a local maximum in the rate of 
increase of vorticity is produced Such maxima were spatially 
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correlated with regions of subsequent severe storm development The 
contours of relative vorticity served to define the raesoscale nature of 
the dominant forcing mechanism (fronts, dry line, meso-low) The 
tomadic Neosho storm developed in an area of local maximum of relative 
vorticity The sensitivity of the vorticity field to perturbations in 
the density and number of wind vectors illustrated an important limita- 
tion in tracking winds associated with real weather systems This 
sensitivity would also seem to indicate the interdependence of the 
dynamic and moisture fields with the type of objective analysis scheme 
utilized Future studies should include a stnicture function analysis 
of the wind fields to determine the "optimal" analysis scheme that 
should be used to minimize the noise and error. 

A second technique of observation was applied after the initiation 

of a cloud element In the early phase of cloud growth (first 30 

minutes) the temporal change of cloud area was monitored. Infrared 

cloud radiances provided fields of equivalent blackbody temperature. 

Areal changes of representative isotherms are then proportional to 

cloud-top divergence, and thus indirectly to the vertical motion 

occurring beneath the cirrus anvil. Severe storms were found to have 

-3 -1 

initial growth rates of 2 x 10 s or greater, while non-severe storms 
grew out at a significantly slower rate. Neither the type of severe 
occurrence nor the time of occurrence on the ground could have been 
reliably predicted based on the cloud-top infrared signatures. The 
small sample size, however, limits the generalization of results. 
Relative cloud-top temperatures were another good indicator of storm 
severity, with clouds 10**C or colder than their neighbors generally 
becoming severe. Minimum temperature changes of 1.5'’C min ^ or greater 
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were also indicative of storm severity Some limitations in measuring 
growth rates in the infrared included the coarse resolution and finite 
response time of the sensor. The ne.ed for multispectral imaging of 
severe storms seems apparent. Such a multispectral technique might 
consist of cloud-top temperatures from infrared radiances, visible 
channel radiances to detect overshooting domes and a polarized micro- 
wave sensor to detect precipitation The need for 5-minute interval 
coverage is also of importance because conventional 30-minute imagery 
cannot resolve the mesoscale and convective scale processes necessary 
to produce a short-term "point" forecast 

A detailed comparison of radar and satellite data was possible for 
the Wewoka hailstorm Excellent correlation existed between temporal 
variation of total echo area and blackbody temperature contours For 
intense echo area C> 30 dbz) no corresponding temperature contour could 
be identified. Using hourly NWS radar summaries, results indicated that 
for storms penetrating the tropopause, large scatter existed between 
radar-observed tops and those computed from blackbody temperatures and 
an environmental sounding. The problem is likely one of resolution, 
the infrared sensor cannot detect small, high-rising turrets nor small, 
intense echo areas beneath the cirrus anvil 

The utilization of a computer- linked, interactive video system to 
process and display digital imagery must also be an integral part of 
any "real-time" application of these techniques The computer supplies 
the rapid processing and infinite patience necessary in displaying 
satellite data, while the satellite meteorologist can apply the known 
principles of atmospheric motion to produce as an integrated result a 
better severe local storm forecast It is in this way that the real 
benefits of quantitative satellite measurements can be realized 
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APPENDIX A 

ALL DIGITAL VIDEO IMAGING SYSTEM FOR ATMOSPHERIC RESEARCH 

The All Digital Video Imaging System for Atmospheric Research 
(ADVISAR) IS a sophisticated blend of hardware and software components 
designed to let the user interact subjectively with his data set It 
is particularly suited to the satellite meteorologist who often has a 

g 

temporal sequence of images, each of which may have as many as 10 
"words" or discrete pieces of information. The overall system is 
described in terms of its five component subsystems 

1 Computer hardware and peripherals 

2 Video display hardware and control 

3. Digitizing hardware and control 

4, System executive software 

5 Applications software 

A pictorial view of the ADVISAR is shown in Figure A-1 Individual 

hardware components are lettered and are listed in Table A-1. Acces- 
sible data storage consists of eight 512 x 512 rapid refresh video 
memories. So for instance, eight consecutive images can be fashioned 
into a temporal loop and displayed under computer or manual control 
Dynamic enhancement tables allow the user to select a range of data 
values and to adjust the input/output "lookup" according to his specific 
needs An electronic cursor allows the user access to every element 
in the 512 x 512 array Modular programming permits a relatively naive 
user to write his own meteorological parameter extraction technique 
(i.e , wind tracking, growth rates, bispectral imaging) which can then 
be run in conjunction with existing applications software Such a 
system is limited solely by the imagination of its users. 
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Table A-1. Ail Digital Video Imaging System for 
Atmospheric Research (ADVISAR) 


A DIGITIZING SUB-SYSTEM 

B 108 CONTROLLER AND B/W MONITOR 

C DATA COLOR EDGE ENHANCER 

D SCAN CONVERTER 

E HP 2100-A COMPUTER 

F DIGITAL DISK 

G 9 TRACK DRIVE 

H 7 TRACK DRIVE 

I CARD READER 

J SYSTEM TELEPRINTER AND KEYBOARD 

K JOYSTICK 

L COLOR DISPLAY MONITOR 

M OPERATOR CRT AND KEYBOARD 
N SOLID STATE RAPID VIDEO REFRESH MEMORIES 
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APPENDIX B 
GLOSSARY OF TERMS* 


conditional inst-abi-lity — the state of a column of air in the atmosphere 
when Its lapse rate of ten^ierature is less than the dry-adiabatic 
lapse rate but greater than the saturation-adiabatic lapse rate 

convection — atmospheric motions that are predominantly vertical, 
resulting in vertical transport and mixing of atmospheric 
properties 

difluence — the rate at which adjacent flow is diverging along an axis 
oriented normal to the flow at the point in question 

entrainment — the mixing of environmental air into a pre-existing 

organized air current so that the environmental air becomes part 
of the current. 

frontogenesis — an increase in the horizontal gradient of an air mass 

property, and the development of the accompanying features of the 
wind field that typify a front. 

infrared radiation — electromagnetic radiation lying in the wavelength 
interval from about 0 8 microns to an indefinite upper boundary 
usually set at 100 microns (0.01 cm) 

lapse rate — the decrease of an atmospheric variable with height, the 
variable being temperature, unless otherwise specified 

level of free convection (LFC) — the level at which a parcel of air 
lifted dry-adiabatically until saturated and saturation- 
adiabatically thereafter would first become warmer than its 
surroundings in a conditional unstable atmosphere 

lifted condensation level (LCL) — the level at which a parcel of moist 
air lifted diry-adiabatically would become saturated 

lifted index — a measure of the potential instability computed by lifting 
the mean moisture in the lower 3,000 feet of the atmosphere moist- 
adiabatically to 500 millibars and subtracting the temperature at 
this point from the reported 500 millibar free-air ten^erature 
\ 

mesoscale meteorology — this field is concerned with the detection and 
analysis of the state of the atmosphere as it exists between 
meteorological stations 

mixing ratio — the dimensionless ratio of the mass of water vapor to the 
mass of dry air. 


*Definitions from Glossary of Meteorology , edited by Ralph E. Huschke, 
American Meteorological Society, Boston, MA (1959) 638 p 
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precipitable water- -the total atmospheric water vapor contained in a 

vertical column of unit cross-sectional area extending between any 
two specified levels 

shortwave — with regard to atmospheric circulation, a progressive wave in 
the horizontal pattern of air motion with dimensions of 1000 to 
2500 km 

Showalter index — a measure of the local static stability of the 
atmosphere, expressed as a niraierical index 

synoptic scale — the study and analysis of weather information on a scale 
of 1000 to 2500 km 

thickness — the vertical depth, measured in geometric or geopotential 

units, of a layer in the atmosphere bounded by surfaces of constant 
pressure 

tropopause — the boundary between the troposphere and stratosphere, usually 
characterized by an abrupt change of lapse rate 

vapor pressure- -the partial pressure exerted by molecules of water vapor 
in the atmosphere 

voxticity — a vector measure of local rotation in a fluid flow, defined 
mathematically as the curl of the velocity vector 
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